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We present an experimental and theoretical study of a scalar atomic magnetometer using an 
oscillating field-driven Zeeman resonance in a high-density optically-pumped potassium vapor. 
We describe an experimental implementation of an atomic gradiometer with a noise level below 
10 fT Hz" 1 / 2 , fractional field sensitivity below 10 -9 Hz" 1//2 , and an active measurement volume of 
about 1.5 cm 3 . We show that the fundamental field sensitivity of a scalar magnetometer is deter- 
mined by the rate of alkali-metal spin-exchange collisions even though the resonance linewidth can 
be made much smaller than the spin-exchange rate by pumping most atoms into a stretched spin 
state. 
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PACS numbers: 07.55.Ge,32.80 Bx,33.35+r,76.60-k 



I. INTRODUCTION 

High-density hot alkali-metal vapors are used in such 
vital metrology applications as atomic clocks [l[ and 
magnetometers @, H, 0|- In these applications the res- 
olution of frequency measurements of the hyperfine or 
Zeeman resonance can be improved by increasing the 
density of alkali-metal atoms until the resonance be- 
gins to broaden due to alkali-metal spin-exchange (SE) 
collisions. Such broadening can be completely elimi- 
nated for Zeeman resonances near zero magnetic field 
[H B 0] • The broadening of the hyperfine and Zeeman 
resonances at a finite magnetic field can be reduced by 
optically pumping the atoms into a nearly fully polar- 
ized state [8j, y, [lOj . These techniques have been used to 
demonstrate clock resonance narrowing [9( and have led 
to significant improvement in the sensitivity of atomic 
magnetometers [nj and to their application for detec- 
tion of magnetic fields from the brain [l2j and nuclear 
quadrupole resonance signals from explosives [l3l |. How- 
ever, the effects of SE collisions on the fundamental sen- 
sitivity of magnetometers operating in a finite magnetic 
field and on atomic clocks have not been analyzed in 
detail. Here we study experimentally and theoretically 
the effects of SE collisions in an atomic magnetometer 
opera ting in geomagnetic field range. It was shown in 
[g, H, llfj| | that in the limit of weak excitation the Zee- 
man and hyperfine resonance linewidths can be reduced 
from Aw ~ R se , where R se is the alkali-metal SE rate, 
to Aw ~ (RgeRsd) 1 ^ 2 , where R s d is the alkali-metal spin- 
destruction rate, by pumping most of the atoms into the 
stretched spin state with maximum angular momentum. 
Since for alkali-metal atoms R s d <C R se (for example, for 
K atoms R s d ~ 10" 4 i? se ), this technique can reduce the 
resonance linewidth by a factor of 10 — 100. However, 
the frequency measurement sensitivity depends not only 
on the linewidth but also on the amplitude of the spin 
precession signal, and the optimal sensitivity is obtained 
for an excitation amplitude that leads to appreciable rf 
broadening. In this paper, we study rf broadening in the 
presence of non-linear evolution due to SE collisions and 
find that the fundamental limit on sensitivity is deter- 



mined by R se even when most atoms are pumped into 
the stretched spin state and the resonance linewidth is 
much narrower than R se . We derive a simple relation- 
ship for the ultimate sensitivity of a scalar alkali-metal 
magnetometer, which also applies qualitatively to atomic 
clocks. We find that the best field sensitivity that could 
be realized with a scalar alkali-metal magnetometer is 
approximately 0.6 fT/Hz 1 / 2 for a measurement volume 
of 1 cm 3 . 

Scalar magnetometers measure the Zeeman resonance 
frequency proportional to the absolute value of the mag- 
netic field and can operate in Earth's magnetic field. 
They are important in a number of practical applications, 
such as mineral exploration [14|, searches for archeolog- 
ical artifacts [l5| and unexploded ordnance [l6j], as well 
as in fundamental physics experiments, such as searches 
for a CP-violating electric dipole moment 0. Some of 
these applications require magnetometers that can mea- 
sure small (~ fT) changes in geomagnetic-size fields with 
a fractional sensitivity of 10" 10 — 10" 11 . Existing sensi- 
tive scalar magnetometers use large cells filled only with 
alkali- metal vapor and rely on a surface coating to reduce 
relaxation of atoms on the walls 0, d, 0] • Here we use 
helium buffer gas to reduce diffusion of alkali atoms to 
the walls, which also allows independent measurements 
of the magnetic field at several locations in the same 
cell [ll[. We present direct measurements of the mag- 
netic field sensitivity in a gradiometric configuration and 
demonstrate noise level below 10 fT Hz" 1//2 in a 10" 5 T 
static field (1 part in 10 9 ) using an active measurement 
volume V ~ 1.5 cm 3 . A small active volume and the 
absence of delicate surface coatings opens the possibil- 
ity of miniaturization and batch fabrication [T^ | of ultra- 
sensitive magnetometers. The best previously-reported 
direct sensitivity measurement for a scalar magnetome- 
ter, using a comparison of two isotopes of Rb occupying 
the same volume V — 180 cm 3 , had Allan deviation that 
corresponds to sensitivity of 60 fT Hz" 1 / 2 and fractional 
sensitivity of 5 x 10" 8 Hz" 1 / 2 [lj|. Theoretical estimates 
of scalar magnetometer sensitivity based on photon shot 
noise level on the order of 1 fT Hz" 1 / 2 have been reported 
in cells with V ~ 1000 cm 3 ||, [|. 



2 



We rely on a simple magnetometer arrangement using 
optical pumping with circularly-polarized light parallel 
to the static magnetic field B z , excitation of spin coher- 
ence with an oscillating transverse magnetic field 2Bi, 
and detection of spin coherence by optical rotation of a 
probe beam orthogonal to the static field. RF broad- 
ening of magnetic resonance is usually described by the 
Bloch equations with phenomenological relaxation times 
Ti and T 2 19]. Since SE collisions generally cause non- 
linear spin evolution, such a description only works for 
small spin polarization 20]. To study the general case 
of large polarization and large rf broadening we per- 
formed measurements of resonance lineshapes in K va- 
por for a large range of SE rates, optical pumping rates, 
and rf excitation amplitudes. We also developed a pro- 
gram for numerical density matrix modeling of the sys- 
tem. To understand the fundamental limits of the mag- 
netometer sensitivity, we derive an analytical result that 
gives an accurate description of magnetometer behavior 
in the regime R se 3> R op ^> R s d, where R op is the opti- 
cal pumping rate, applicable to high density alkali-metal 
magnetometers with high spin polarization. In the limit 
of high polarization, we find an implicit equation for the 
transverse spin relaxation Ti that can be solved to cal- 
culate polarization P as a function of rf field detuning 
and other parameters. In this limit, the system is well- 
described by the solutions to the familiar Bloch equa- 
tions, with T2 varying as a function of polarization and 
rf field de-tuning. This modified Bloch equation model 
reproduces the non-Lorentzian resonance lineshape from 
the full density matrix simulation and the experimental 
rf broadening data and allows us to set analytical lim- 
its on the magnetometer sensitivity. The same approach 
can also be easily applied to other alkali metal atoms 
with different nuclear spin values and to hyperfine clock 
transitions. 

This paper is organized as follows: Sec.[TT]describes the 
experimental setup and presents measurements of mag- 
netic field sensitivity and other experimental parameters. 
Sec. IIIII presents a theoretical description of the magne- 
tometer signals. Sec. IIVI gives expressions for the fun- 
damental sensitivity of the magnetometer and compares 
this theoretical result to our high-sensitivity magnetome- 
ter measurements. 



II. EXPERIMENTAL MEASUREMENTS 

A. Measurement apparatus 

The scalar magnetometer, diagrammed in Fig. [1] is 
built around a Pyrex cell containing potassium in natural 
abundance, 2.5 amg of 4 He to slow atomic diffusion, and 
60 Torr of N2 for quenching. For characterization, the cell 
was heated to varying temperatures using a hot air oven. 
For the most sensitive magnetometry measurements, the 
cell was heated with pairs of ohmic heaters (wire meander 
in Kapton sheet) oriented to cancel stray fields and driven 




FIG. 1: Schematic of the experimental apparatus. The cell 
(3x3x4 cm 3 with the larger dimension perpendicular to 
lasers) is heated inside a boron nitride oven and placed in a 
glass vacuum enclosure pumped out to 0.5 Torr. Coils inside 
6-layer magnetic shields allow application of magnetic fields 
and gradients. The gradiometer measurement is obtained by 
imaging the probe beam onto two-element photodiodes. The 
signals of the two balanced polarimeters are subtracted at the 
lock-in. 



at 27 kHz. A circularly polarized pump beam at the Dl 
resonance polarizes the K atoms along the z-direction. 
The x component of atomic spin polarization is mea- 
sured using optical rotation of a linearly-polarized beam 
as determined by a balanced polarimeter. Two-segment 
photodiodes were used on each arm of the polarimeter 
to make a gradiometer measurement. A constant bias 
field B z is applied parallel to the pump laser. An os- 
cillating rf field 2Bi is applied in the y direction with 
its frequency tuned to the Zeeman resonance given by 
cj = jB z = g s [i B B z /(2I+l)h = 27rx(700kHz/G) J B z for 
potassium atoms. The polarimeter measurement is read 
through a lock-in amplifier, tuned to the rf frequency. 
The lock-in phase is adjusted to separate the resonance 
signal into symmetric (in-phase) absorption and antisym- 
metric (out-of-phase) dispersion components. Exactly on 
resonance the dispersive part of the signal crosses zero. 
The magnitude of the local field is determined by the fre- 
quency of this zero-crossing and changes in the dc mag- 
netic field are registered as deviations from zero of the 
dispersive signal. 



B. Noise measurements with a high sensitivity 
atomic magnetometer 

Magnetometer noise is read on the dispersive compo- 
nent of the lock-in reading. The conversion of the voltage 
noise to magnetic field noise depends on the slope as a 
function of the magnetic field or frequency of the disper- 
sion curve. The tunable parameters of the experiment 
were adjusted to maximize the dispersion curve slope. 
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FIG. 2: Absorptive (open symbols) and dispersive (closed 
symbols) components of the magnetic resonance polarization 
rotation signal at 1 /xT (squares), 10 /xT (triangles), and 26 piT 
(circles). Solid lines show Lorentzian fits to the data. These 
data were recorded at the same time and under the same 
experimental conditions as the high sensitivity magnetometer 
measurements. 



The pump beam (20-40 mW) was imaged on an area of 
roughly 3 x 1.5 cm 2 across the cell. A probe beam cross 
section of 1.2 x 1.2 cm 2 was defined by a mask with total 
power of 10 mW and the wavelength detuned by about 
100 GHz from the Dl resonance. After passing through 
the cell and the polarizing beam splitter the probe beam 
was imaged onto two-segment photodiodes. For the most 
sensitive measurements, the amplitude of the oscillating 
rf field was about 19 nT. Magnetic field sensitivity was 
measured for three values of B z : 1 /xT, 10 /xT, and 26 
/xT. The cell was heated to approximately 150°C, yield- 
ing an atomic density of n = 6.4 x 10 12 cm~ 3 . The po- 
larimeter signals were measured with a lock-in amplifier 
(Stanford Research Systems SR830 for 1 /iT and 10 /xT 
measurements, SR844 for the 26 /xT measurement). The 
lock-in internal reference generated the rf field and the 
time constant was set to 100 (is. The resonance line- 
shapes obtained by varying the rf frequency are shown 
in Fig. [21 The pump power and rf amplitude are ad- 
justed to optimize the slope of the dispersion signal for 
a given probe beam power. At the parameters that opti- 
mized the magnetometer sensitivity, the resonance curves 
are well-described by Lorentzian lineshapes with similar 
half-width at half maximum (HWHM) for absorptive and 
dispersive components of ~ 220 Hz for 1 /iT and 10 /xT 
and 265 Hz for 26 /xT. The amplitude and width of the 
optical rotation signal was found to be nearly indepen- 
dent of the static magnetic field values over the range of 
our measurements. The field B z was generated using a 
custom current source, based on a mercury battery volt- 
age reference and a FET input stage followed bya con- 
ventional op-amp or a transistor output stage [2jJ . The 
fractional current noise was less than 2 x 10" 8 Hz 1 / 2 at 10 
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FIG. 3: Noise spectra for 1 /xT (a) and 10 /xT (b). Shown are 
single channel spectra (black line with crosses), two-channel 
difference (gradiometer) spectra (black solid line), and the 
measured electronic and optical noise (gray solid line) ob- 
tained by blocking the pump beam. The dashed black line 
marks the 14 fT/Hz 1 ^ 2 level. Magnetic field noise increases 
at higher frequencies due to correction for the finite band- 
width of the magnetometer. 



Hz, about 10 times better than from a Thorlabs LDC201 
ULN current source. Low-frequency (< 10 Hz) optical 
rotation noise was reduced by an order of magnitude by 
covering the optics with boxes to reduce air convection 
that causes beam steering. The oven and laser beams 
within the magnetic shields were enclosed in a glass vac- 
uum chamber to eliminate air currents. 

Probe beam position was adjusted to equalize the pho- 
todiode signals for the two polarimeters within 2%. The 
gradiometer measurements reduced by more than an or- 
der of magnitude the noise from the B z current source 
as well as pump intensity and light shift noise. By ap- 
plying a calibrated magnetic field gradient, we found the 
effective distance between the gradiometer channels to 
be ~ 3.5 mm, much larger than the K diffusion length 
in one relaxation time (DT2) 1 / 2 ~ 0.1 mm, so the two 
measurements are independent. 

The magnetic field data were acquired from the dis- 
persive lock-in signal for 100 sec with a sampling rate of 
2 kHz. The FFT of the data was converted to a magnetic 
field noise spectrum using a frequency calibration of the 
dispersion slope and corrected for the finite bandwidth 
of the magnetometer. The bandwidth was found to be 
close to the Lorentzian HWHM for all values of B z . The 
magnetic noise spectra at 1 /xT and 10 /xT are shown 
in Fig. [3J At 1 /xT, single channel measurements were 
limited by lock-in phase noise, while at 10 /xT they were 
limited by current source noise. The noise in the differ- 
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ence of the two channels was limited almost entirely by 
photon shot noise at higher frequencies and reached be- 
low 14 fT/Hz 1/2 , corresponding to less than 10 fT/Hz 1/2 
for each individual magnetometer channels. With the 
pump beam blocked, the optical rotation noise reached 
the photon shot noise level. Low frequency noise was 
most likely due to remaining effects of convection. At 
26 juT, the gradiometer had a sensitivity of 29 fT/Hz 1//2 , 
limited by lock-in phase noise and imperfect balance be- 
tween gradiometer channels. 



C. Magnetic resonance measurements 

To analyze the magnetometer behavior and predict 
the theoretical sensitivity of the device, we focus on the 
shape of the magnetic resonance curves. The basic be- 
havior of the resonance signals can be understood us- 
ing phenomenological Bloch equations (BE), which pre- 
dict a Lorentzian resonance lineshape. Though the BE 
cannot describe the whole physics in the case of rapid 
spin-exchange collisions, they do provide a convenient 
phenomenological framework for qualitative understand- 
ing of the resonance lineshape including the effects of rf 
broadening. Using the rotating wave approximation, the 
solution of the BE (see, for example [lj|) in a frame ro- 
tating about the z-axis is: 



P = 



Pz = 



1 + {AluT 2 ) 2 + ( 7 Bi) 2 TiT 2 

jBiT 2 

l + (AwT 2 ) 2 + ( 7 S 1 ) 2 r 1 T 2 

1 + {AuT 2 f 
1 + (AwT 2 ) 2 + (7Bi) 2 TiT 2 ' 



Pa 



(1) 
(2) 
(3) 



Here we introduce the in-phase Pj and out-of-phase Pi 
components of the transverse polarization in the rotat- 
ing frame and the longitudinal polarization P z . In the lab 
frame, we measure P x = PjCos(ojt) + PiSin(ujt), and we 
tune the lock-in phase to separate the absorptive Pj from 
the dispersive Pi. T\ and T 2 are constant phenomenolog- 
ical relaxation times, Pq is the equilibrium polarization, 
B\ is the amplitude of the excitation field in the rotating 
frame, given by B y — 2_B 1 cos(tji) in the lab frame. The 
detuning Auj = u> — ujq is the difference between the rf 
frequency u> and the resonant frequency loq, which is the 
Larmor frequency in the applied dc field B z . The depen- 
dencies of Pi and Pj on frequency are Lorentzian, with 
the HWHM 
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The increase in the width due to the presence of exci- 
tation field Bi is the basic phenomenon of rf resonance 
broadening. The slope, at resonance, of the dispersive 
component of the signal dPi/duj(Auj = 0) is given by 



(IlU 



1 B 1 T 2 2 



The slope has a maximum at an excitation field B\ 
1/(7V^2)= 



P T 1 2 3/2 /(2T 1 I 



1/2, 



(6) 



The accuracy of the simple Bloch equation theory de- 
pends on the contribution of spin-exchange relaxation to 
the linewidth. If the temperature is low, then the broad- 
ening due to optical pumping can exceed SE broadening, 
and the Bloch equation theory will be quite accurate. 
Additionally, if spin polarization is low, SE broadening 
will not depend significantly on the polarization and the 
excitation field, so the transverse relaxation time T 2 will 
be almost constant; in this case the BE solution is also 
valid. However, we are primarily interested in the regime 
of high spin-exchange rate and high spin polarization, 
where the magnetometer is most sensitive. 

To understand the effects of SE broadening, we com- 
pared the lineshape predicted from the BE to the mea- 
sured resonance lineshape of the magnetometer signal at 
a frequency of 80 kHz. We recorded the magnetic reso- 
nance curves for different values of rf excitation ampli- 
tude, pump laser intensity, and cell temperature. We find 
that the lineshape of the resonance remains reasonably 
close to a Lorentzian and the in- and out-of-phase lock-in 
data from resonance measurements were fit to the absorp- 
tive and dispersive Lorentzian profiles, allowing for some 
mis-tuning of the lockin phase. It can be seen from the 
BE that relative amplitudes of absorptive and dispersive 
components can differ substantially from those expected 
from a complex Lorentzian l/[i(u> — u>o)+T] in the regime 
of large rf broadening. Moreover, due to SE effects, the 
absorption and dispersion widths for the same experi- 
mental conditions can also differ. Thus, a total of five 
parameters were used for each resonant curve: the reso- 
nant frequency, and the respective amplitudes and widths 
of the absorptive and dispersive signal components. An 
example of the results for the resonance linewidths as a 
function of the magnitude of the rf field is shown in Fig.|4j 
It can be seen that rf broadening is greater than what is 
expected from the BE. Moreover, the absorptive and dis- 
persive parts of the resonance have different widths as 
the rf amplitude in increased. These are signatures of 
the SE broadening that require modifications of the BE 
description. 

In the regime of small rf broadening we verified that 
the absolute size of the lock-in signal is in agreement with 
Bloch equations. The optical rotation signal detected by 
the lock-in is given by 



lr e cfnD{ V )P z 
2^2 



(7) 



1 + { 7 B 1 ) 2 T 1 T 2 



(5) 



where D(v) = (v — vq)/{{v — vo) 2 + T 2 ) is the optical 
dispersion profile of the Dl resonance line with linewidth 
T v and oscillator strength /, n is the density of atoms, 
I is the length of the cell in the direction of the probe 
beam, and r e = 2.8 x 10 -13 cm is the classical electron 
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FIG. 4: The linewidths of Lorentzian fits to the experimen- 
tal data for absorption (solid points) and dispersion (open 
points) components of the magnetic resonance in K vapor at 
140°C. The dotted line is the prediction for rf broadening of 
the linewidth from Bloch equations with constant T\ and T2, 
the solid and dashed lines are results of Lorentzian fits to 
absorption and dispersion lineshapes obtained using modified 
BE with variable T2 discussed discussed in the text. Here 
R se = 5100 s" 1 and R ad = 24 s _1 are fixed from indepen- 
dent measurements, while R op — 840 s _1 is adjusted to fit 
the measured linewidth at low rf amplitude. 



radius. Here we take into account the fact that lock-in 
output measures the r.m.s. of an oscillating signal. The 
length I is determined by the dimensions of polarized 
vapor illuminated with the pump beam. Near the edges 
of the cell the pump beam is distorted, reducing I below 
the inner dimensions of the cell. We varied the width 
of the pump beam to find that the largest pump width 
for which the signal still increases is about 2 cm. For 
this value of / the absolute signal size was in agreement 
with Bloch equations to within 15%. So the volume of the 
polarized atomic vapor participating in the measurement 
is about 2 x 1.2 x 0.6 cm 3 ~ 1.5 cm 3 . 



D. Relaxation rates 



TABLE I: Comparison of the measured SD rates R™ d and 
SE rates R^l from fits of the resonance linewidth at low field 
with corresponding rates R^ and R^ 1 calculated from colli- 
sion cross sections and the density of K metal determined by 
optical absorption. 
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140 


22±2 
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3.8 
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43±6 


23±2 


10±0.5 


8.1±0.2 


6.7 


160 


31±3 


25±2 


14±1.0 


13.9±0.2 


11.4 



was fit to a sum of two Lorentzians taking into account 
the counter-rotating component of pump rate modulation 
[22| . At low magnetic field, when the Zeeman frequency 
is much smaller than the SE rate, SE broadening depends 
quadratically on the magnetic field. The spin-destruction 
rate is obtained from extrapolation of the width to the 
zero-field limit. From these fits of the resonant frequency 
and linewidth we determined the spin-exchange rate R se 
and the spin-destruction rate (SD) R s d, which are listed 
in Table U for the same cell at several temperatures. The 
error bars are estimated from fits to different sets of the 
data. In addition, we determined the density of K atoms 
by scanning the DFB probe laser across the the optical 
absorption profile of the Dl resonance. The alkali densi- 
ties calculated from the integral of the absorption cross- 
section using known oscillator strength (/ = 0.34) and 
cell length (1 = 3 cm) are also shown in Table Q] The den- 
sity is approximately a factor of 2 lower than the density 
of saturated K vapor at the corresponding temperature, 
as we find is common in Pyrex cells, probably due to slow 
reaction with glass walls. The alkali-metal SE rate can be 
calculated from the measured density using known K-K 
spin-exchange cross-section ose — 1-78 x 10~ 14 cm 2 [23| 
and is in good agreement with direct measurements. The 
spin destruction rate R s d can also be calculated using 
previously measured spin-destruction cross-sections for 
K-K, K-He and K-N2 collisions [24| and gas composi- 
tion in the cell (2.5 atm of 4 He and 60 torr of N2). We 
also include relaxation due to diffusion to cell walls. Er- 
rors on the rates calculated from the densities are esti- 
mated from uncertainty in the cross-section and in the 
gas pressures in the cell. Our direct measurements of 
the spin-destruction rate are reasonably consistent with 
these calculations. 



A number of independent auxiliary measurements were 
preformed to find the relaxation rates of the alkali-metal 
spins to be used for detailed modeling of SE effects. Spin 
exchange and spin destruction rates can be determined 
by measuring the width of the Zeeman resonance in a 
very low field using low pump and probe laser inten- 
sity For these measurements, the magnetic field was 
perpendicular to the plane of the lasers and the pump 
laser intensity was modulated near the Zeeman reso- 
nance. The signal as a function of modulation frequency 



III. MODEL OF MAGNETOMETER 
DYNAMICS 

We first model the dynamics of the system using nu- 
merical evolution of the density matrix to accurately de- 
scribe the effects of SE relaxation. To provide more quali- 
tative insight and estimate the fundamental limits of sen- 
sitivity we also develop a semi-analytical description, a 
modification of the BE, that provides a good approxi- 
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mation to the numerical solutions in the regime of high 
spin-exchange rate. 



A. Density matrix equations 

The spin evolution can be accurately described by the 
solution of the Liouville equation for the density matrix. 
The time-evolution of the density matrix p(t) includes 
hyperfine interaction, static and rf field interactions, op- 
tical pumping, spin relaxation processes and non-linear 
evolution due to alkali- metal spin-exchange collisions. In 
the presence of high density buffer gas when the ground 
and excited state hyperfine structure of the alkali-metal 
atoms is not resolved optically, the density matrix evolu- 
tion is given by the following terms [251 ]: 



[B-S,p} + — — (8) 



T, 



sd 



<p(l + 4(S)-S)-p 



R op [<p{l + 2s- S)- p] 



Here, A^f is the hyperfine coupling, I is the nuclear spin 
and S is the electron spin operator. The Bohr magne- 
ton is [Ib and g s is the electron g-factor, B is the exter- 
nal magnetic field including static and oscillating compo- 
nents, ip is the purely nuclear part of the density matrix 
[25l | . and s is the spin polarization of the pump beam. 
We evaluate the density matrix in the \F,m) basis and 
focus on the regime of relatively low static magnetic field, 
where the non-linear Zeeman splitting given by the Breit- 
Rabi equation is small. To simplify numerical solution of 
the non-linear differential equations we neglect hyperfine 
coherences and make the rotating wave approximation 
for Zeeman spin precession, 

(F,m\p(t)\F',m') = 5 F F , (F, m\p'(t)\F, m')e zuj{m ' ~ m)t . 

(9) 

Here u> is the frequency of rf excitation field tuned near 
the Zeeman resonance and (F,m\p'(t)\F,m') is the den- 
sity matrix element in the rotating frame, evolving on 
a time scale on the order of spin relaxation rates that 
are much slower than the Zeeman spin precession fre- 
quency. With this approximation it is necessary to con- 
sider only 21 elements of the density matrix for / = 3/2 
using the symmetry of the off-diagonal components. In 
the rotating frame without loss of generality we pa- 
rameterize the density matrix as p'(t) = p' ST (P,9,<j)) + 
p[. Here p' ST (/3,9,<j>) is a spin-temperature distribution 
(F,m\p' ST (j3)\F, m') cx e^ m that is rotated by an angle 
9 from the z axis into the j direction of the rotating 
frame and an angle <fi around the z axis. p\ is a den- 
sity matrix describing deviations from spin-temperature 
distribution, which are small because the spin-exchange 
rate is much larger than all other rates. For a given 
value of the spin temperature j3 and angles 9 and (j> the 
expectation value of (S) = Tr[Sp' ST ((3, 6, 0)] is used in 
the spin-exchange term of the density matrix evolution 



equations, reducing them to a set of linear first order dif- 
ferential equations for the perturbation matrix p' x . The 
steady-state solution for p[ is obtained symbolically in 
Mathcmatica. To obtain a self-consistent solution, f3, 9 
and 4> are adjusted until the steady-state solution for p[ 
satisfies TrfS^] = 0. The self-consistency iteration is 
performed numerically for various values of the optical 
pumping rate and the rf excitation strength and detun- 
ing. 



B. Modified BE 

Though the numerical solutions to the density matrix 
equations give an accurate treatment of the spin dynam- 
ics, it is convenient to develop an analytical model that 
can describe the asymptotic behavior of the system in the 
regime of high spin-exchange rate. Here we focus on the 
regime of light-narrowing [aEIJ, with R se 3> R op ^> R s d, 
which also implies that P is close to unity. For weak rf 
excitation an analytic expression for T2 under these con- 
ditions has been obtained in (§, LLO, k3 , 



J_ _ Rop Rse n _ p v 

T 2 ~ 4 + 5 1 z) - 



(10) 



The coefficients in this expression depend on the nuclear 
spin / and on the size of the nonlinear Zeeman splitting 
relative to the spin-exchange rate [lOj • Eq. (fTT))) describes 
the case of I = 3/2 and large spin-exchange rate relative 
to the non-linear Zeeman splitting, so all Zeeman reso- 
nances overlap. It is clear from this equation that spin- 
exchange relaxation can be suppressed by maintaining P z 
close to unity. 

To extend this solution to arbitrary rf excitation we 
observe that the relaxation due to spin-exchange and op- 
tical pumping is independent of the direction of spin po- 
larization. Therefore, we can apply Eq. (|10[) in a rotating 
frame with z' axis tilted by an angle 9 from the lab z axis 
and rotating together with P in the presence of a large 
rf excitation field. In doing so we introduce an error due 
to inaccurate treatment of transverse spin components in 
the F = 1 state. Spin precession in F = 1 state occurs 
in the direction opposite to the precession in F = 2 state 
and hence will not be stationary in the rotating frame. 
However, this error is small in the light narrowing regime 
because of two small factors: a) the population in F = 1 
state is small since P is close to unity and most atoms are 
pumped into the stretched state with F = 2 and b)^<Cl 
for rf fields that provide optimal sensitivity to maintain 
P close to unity and hence the transverse components of 
spin are small. 

Using this approximation we then solve BE (Eq. fl|[3])) 
in combination with an equation for Ti as a function of 
polarization 



1 



R 



op 

T 



R, 



■[l - (Pi + Pf + Pt) 



2U/21 



(11) 
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FIG. 5: Comparison of transverse polarization components 
(Pi, Pj) using full numerical density matrix evolution (solid 
points - absorption, open points - dispersion) and modified 
BE (solid line - absorption, dashed line - dispersion) for 
Rse/Rsd = 10 4 , Ro P /Rsd = 200, ~fB 1 /R sd = 100. Lorentzian 
lineshapes are shown with dotted lines for comparison. 



The longitudinal spin-relaxation time is not affected by 
spin exchange and is given by Ti = 4/ (R op + R s d) in the 
limit of high spin polarization [a, [l(| . The equilibrium 
spin polarization in the absence of rf excitation is equal to 
Pq = R op I (R op + R s d)- The resulting algebraic equations 
can be easily solved for arbitrary parameters. However, 
the solution is only expected to be accurate when P re- 
mains close to unity. In Fig. [5] we compare the resonance 
lineshapes obtained with a full numerical density matrix 
and the analytical calculation using modified BE. It can 
be seen that for this case which is well into the asymp- 
totic regime R se R op 3> R s d the analytical results 
agree very well with exact calculations. The lineshapes 
are significantly different from a simple Lorentzian. 



C. Comparison of experimental measurements 
with theory 

The results of the simple BE, numerical calculations 
with the full density matrix equations, and analytical re- 
sults from the modified BE were compared to a large set 
of measurements in various parts of the parameter space. 
One such comparison is shown in Fig. 3J The experi- 
mental data compare well to the BE when a variable T2 
(from Eq. ITT]) is used. Note that not only is the mea- 
sured width greater than that predicted by the simple 
BE (with constant T2) but also, as correctly predicted 
from the analytical theory, the half-width of the absorp- 
tion curve differs from the half-width of the dispersion 
curve. At higher excitation amplitudes, even the modi- 
fied Bloch analysis begins to deviate from the measured 
half widths because the polarization begins to drop. The 
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FIG. 6: Slope of the dispersive part of the resonance curve, 
given as polarization per angular frequency dPi/du. Ex- 
periment: Open circles - R op = 220 s _1 , open triangles - 
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absorbtivity of the vapor also changes as a function of 
the rf excitation and the pumping rate at the location of 
the probe beam is not a constant. This can be taken into 
account by considering the propagation of the pumping 
light through the polarized vapor. Though the width of 
the resonance is a good metric for comparing experiment 
to theory, it is the slope of the dispersive component at 
resonance that is most important for the magnetome- 
ter sensitivity. In Fig. [6] is shown a comparison of the 
measured slopes (from the same data as Fig. 2]) to those 
predicted by the analytical theory. The agreement is gen- 
erally satisfactory. In Fig. [7] we show one example of a fit 
of the measured resonance profile to that predicted from 
the modified BE. As these data show, in the parameter- 
space of interest, the modified Bloch analysis provides a 
good description of the slope and the width of the reso- 
nance, as a function of R op , B\, R s d, and R se . Thus we 
can use these equations to determine the best-achievable 
sensitivity of the scalar magnetometer. 



IV. MAGNETOMETER SENSITIVITY 

For a given slope of the dispersion curve, unavoidable 
noise sources in the system determine the fundamental 
sensitivity of the scalar magnetometer. The calculation 
of the sensitivity follows closely that for an rf atomic 
magnetometer, derived in [Toj ] . For a given polarization 
noise SPi, the resulting field noise SB is 



SB 



SP 



-y\dPi/du\ 



(12) 



There are many sources of technical noise which con- 
tribute either directly to the scalar magnetometer noise 
as in the case of low-frequency magnetic field noise from 
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FIG. 7: Comparison of experimental and theoretical reso- 
nance lineshapes. T = 140 °C, R sd = 24 s" 1 , R se = 5100 



R op — 625 s 



Bi = 20 nT. The model also includes 
a correction for the polarization-dependent absorption of the 
pump beam. 



the current source, or indirectly as in the cases of voltage 
noise of an amplifier, magnetic field noise at high fre- 
quency, vibrations of the optical detection system, and 
pump laser noise. Technical noise can be removed in prin- 
ciple, so it is important to understand the fundamental 
limits that determine the best achievable sensitivity. 



B. Light-shift noise 

The ac Stark shift (light shift) is induced by the probe 
beam, which is tuned off-resonance from the atomic tran- 
sition, when it has a non-zero circular polarization. If 
the probe laser detuning is much larger than the hyper- 
fine splitting, the action of the light on atomic spins is 
equivalent to the action of a magnetic field parallel to the 
light propagation direction. This light-shift field is given 
by (see Eq. (9) of Ref. (3), 



B 



LS 



„cf§ pr s x D(v) 
(21 + 1)jA ' 



(16) 



where s x is the degree of circular polarization of the probe 
beam. Light-shift noise can occur as a result of fluctua- 
tions of intensity, wavelength, or s x . If the probe beam 
is perfectly linearly polarized, fluctuations of the circular 
polarization are due to quantum fluctuations resulting in 
an imbalance between the number of left and right circu- 
larly polarized photons in the probe beam. The spectral 
density of the probe beam spin polarization noise is given 
by 5s x = v/2/$j, r . Substituting this value of s x and ex- 
cluding <& pr by using the pumping rate of the probe beam 
R pr in the limit (y — vq) 3> T nu we get: 



SB 



LS 



yj 2R pr N a b 



(17) 



A. Photon shot noise 

In a balanced polarimeter the polarization rotation 
noise per unit bandwidth due to quantum fluctuations of 
the number of photons received by photodetectors with 
quantum efficiency r\ is given by 



1/V2$ pr 77, 



(13) 



where $ pr is the number of photons per second in the 
probe beam. The noise has a flat frequency spectrum 
and 8(f) is measured in units of rad/Hz 1 / 2 . The same 
level of noise per unit bandwidth will be measured in 
each phase of a lock-in amplifier calibrated to measure 
the r.m.s. of an oscillating signal. The optical rotation 
measured by the lock-in amplifier is given by Eq. ([7]). 

It is convenient to express the photon flux <f> pr in terms 
of the pumping rate of the probe beam R pr : 



r e cf§ pr {T v /A) 



(14) 



where A is the cross-sectional area of the probe beam. If 
the probe laser is detuned far from resonance, \v — vo\ S> 
r„, then one can express photon- atom interactions in 
terms of the number of absorption lengths on resonance 
N a b = r e cfnl/T v . Using Eq. (fT2"]) and magnetometer 
volume V = IA we find the magnetic field noise due to 
photon shot noise is given by 

SBp S = 2 (15) 

^dPi/ckjyNabRprnVr] 



This effective field noise {B x s <C Si) causes polarization 
noise by rotating the P z component into the direction 
of the primary signal Pj. The amount of polarization 
noise in Pj induced by the light-shift field is proportional 
to the spin coherence time T 2 . Using simulations of BE 
with noise terms one can verify that 



Pt S = jB^ s T 2 P z /x/2 



(18) 



where a factor 1 / y2 appears because only the component 
of the light-shift field that is co-rotating with the spins 
contributes to the noise. We get the following contribu- 
tion of the light shift to the noise of the magnetometer 



SB 



LS 



PzVR pr N ab T 2 " 



(19) 
~ 1. 



4jVnV\dPi/doj\ 
In most cases of interest here one can assume that P z 

C. Spin projection noise 



The spin-projection noise occurs as a result of quantum 
fluctuations in the components of atomic angular momen- 
tum. We consider the case when the polarization is close 
to unity and most atoms are in F = 2 state. Using the 
fundamental uncertainty relationship SF x SF y > hF z /2 
with (F z ) ~ 2 one can show [l(| that the polarization 
noise per unit bandwidth is given by 



SPi = xjT 2 /N, 



(20) 
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where N is the total number of atoms. The spin projec- 
tion noise depends only weakly on absolute spin polariza- 
tion; for K atoms with I = 3/2, it increases by y/3/2 for 
unpolarized atoms. The resulting magnetic field noise in 
the scalar magnetometer is given by 



SB 



SP 



f\dPi/du\ 



(21) 



D. Optimization of fundamental sensitivity 

Combining all the noise contributions we obtain the 
following equation for the magnetometer sensitivity 



5B 



L dui J 



T^RprNab 

16 



RprNabV' 



(22) 



The first term describes spin projection noise, the second, 
the light shift of the probe beam, and the third, photon 
shot noise. 

To find the fundamental limit of the sensitivity we as- 
sume that Nab can be adjusted separately, for example by 
increasing the length of the sensing region in the probe 
direction while keeping the volume constant, or changing 
the buffer gas pressure. We find that the optimal optical 
length is equal to N a b — 8/ (y/rjT2R pr )- It is always bene- 
ficial to reduce R pr and increase N a b, which will result in 
longer T\ and T% until R pr <C R s d- Under optimal prob- 
ing conditions the fundamental magnetometer sensitivity 
reduces to 



5B 



L du J 



■^(1+7,-1/2). 



(23) 



The best sensitivity is obtained by maximizing 
dPij dui I \fT^. For a given R se and R s d we vary R op 
and B\ and calculate dPi/dto and T 2 using modified BE 
with variable T 2 given by Eq. (fTTj) . We find that for 
R se 3> R s d the maximum value of dPi/cLj/y/T^ is given 

1/2 

by dPi/dw/v^lmax = kR S e , where k = 1.3. This re- 
sult is also verified with the full numerical density matrix 
model. With R se — nva sei the optimal sensitivity of a 
scalar alkali-metal magnetometer is given by 



SB m \ n = 



0.77 va S e(l + y- 1 / 2 ) 
7 V V 



(24) 



Hence we find that for a scalar magnetometer the 
fundamental sensitivity is limited by the rate of spin- 
exchange collisions even though the resonance linewidth 
can be much smaller than the spin exchange rate. Nu- 
merically for age = 1-8 x 10~ 14 cm 2 and r\ = 0.8 we 
find that <5i? m ; n = 0.9 fT/Hz 1 / 2 for an active volume of 
1 cm 3 . Using back-action evasion techniques it is possible 
to make the light shift and photon shot noise contribu- 
tions negligible, but this only improves the sensitivity to 
0.6 fT/Hz 1 / 2 for 1 cm volume. 



If total noise is limited by photon shot noise or by 
technical sources of rotation noise, as was the case in 
our experiment (see Fig. [3]), the sensitivity is opti- 
mized by maximizing the slope on resonance dPi/duj. Us- 
ing the same optimization procedure using modified BE 
and varying R op and B\ one can obtain dPij dio\ msDL — 

\.2R S e^R a d^ 4 ■ ln this case the maximum slope is in- 
creased from R~£ scaling that one would obtain with a 
spin-exchange-broadened resonance from Eq. ©. There- 
fore, light narrowing is useful in reducing the noise in 
scalar magnetometers limited by photon shot noise or 
1// noise, with a maximum sensitivity gain on the order 
of (R Be / Rsd) 1 ^ 4 , which is equal to about 10 for K atoms. 

One can also use Eq. (f2"2"|) to estimate the best sensi- 
tivity possible under our actual experimental conditions. 
In this case the number of absorption length on reso- 
nance N a b is not optimal and the spin relaxation of K 
atoms has additional contribution from collisions with 
buffer gas and diffusion to the walls. For our parame- 
ters corresponding to Fig. [3] (N a b = 2.5, R se = 8700 s _1 , 
R pr ~ 100 s" 1 , Rsd + R P r ~ 130 s _ \ V ~ 1.5 cm 3 and 
7] = 0.24), including losses in collection of probe light 
after the cell), we get optimal sensitivity from Eq. ([22]) 
of 7 fT/Hz 1 / 2 , dominated by photon shot noise. This 
compares well with the measured photon shot noise level 
corresponding to 7 fT/Hz 1 / 2 in each channel. The sen- 
sitivity could be improved by increasing the resonance 
optical depth of the vapor. 



V. CONCLUSION 

In this paper we have systematically analyzed the sen- 
sitivity of a scalar alkali-metal magnetometer operating 
in the regime where the relaxation is dominated by spin- 
exchange collisions. We demonstrated experimentally 
magnetic field sensitivity below 10 fT Hz -1 / 2 with an 
active volume of 1.5 cm 3 , significantly improving on pre- 
vious sensitivities obtained for scalar atomic magnetome- 
ters and opening the possibility for further miniaturiza- 
tion of such sensors. 

We considered the effects of rf broadening in the pres- 
ence of SE relaxation and developed a simple analytic 
model based on Bloch equations with a T 2 time that de- 
pends on rf excitation. The results of the model have 
been validated against a complete numerical density ma- 
trix calculation and experimental measurements. We 
showed that the fundamental sensitivity limit for a scalar 
alkali-metal magnetometer with a 1 cm 3 measurement 
volume is on the order of 0.6-0.9 fT Hz -1 / 2 . In this case 
a reduction of resonance linewidth by optical pumping of 
atoms into a stretched state does not lead to an improve- 
ment of fundamental sensitivity limit. 

It is interesting to compare the scaling of the optimal 
magnetic field sensitivities in various regimes. It was 

shown in [ll[ that near zero field in the SERF regime 

1/2 

the sensitivity scales as a J, , while for an rf magnetome- 
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ter operating in a finite field it scales as {(JseCFsd) 1 ^ [HI- 
In contrast, here we find that the fundamental sensitiv- 
ity limited by spin projection noise for a scalar magne- 
tometer in a finite field scales as ale 2 , i.e. there is no 
significant reduction of SE broadening for optimal con- 
ditions. Since a se is similar for all alkali metals, one can 
expect a similar sensitivity for a Cs or Rb magnetome- 
ter. On the other hand, if one is limited by the photon 
shot noise or technical sources of optical rotation noise, 
which is often the case in practical systems, the mag- 
netometer sensitivity is determined by the slope of the 



dispersion curve. In this case it is improved in the light- 
narrowing regime because the slope of the dispersion res- 
onance scales as <Js e 3 ^ 4 aJ^ 4 , instead of o~ e x for the case 
of spin-exchanged broadened resonance. We expect simi- 
lar relationships, with different numerical factors, to hold 
for atomic clocks operating on the end transitions, since 
To in that case is given by an equation similar to Eq. (|11|) 
[9( . The analytical approach developed in this paper can 
be easily adapted to other alkali atoms by modifying the 
coefficients in Eq. (1111) . This work was supported by an 
ONR MURI grant. 
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